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ABSTRACT 

.> 

The  relation  of  the  total  structure  weight  of  a 
two-dimensional  ring  frame  to  the  stress  distribution  in  each  of 
the  elements  is  investigated.   The  total  weight  of  the  two- 
dimensional  structures  investigated  were  a  minirrun  when  design 
stresses  had  been  attained  in  all  elements  of  the  structure. 

An  expression  is  developed  that  approximates  the  change 
in  stress  distribution  in  a  member  when  additions  or  deletions  cf 
material  are  made  in  the  member  itself  and  its  associated  members. 

This  expression  is  used  as  the  basis  for  a  rational 
design  technique  that  converges  to  a  ring  frame  of  minimum  total 
weight. 
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I.   INTRODUCTION 

The  transverse  structure  of  a  wall  sided  ship  consists 
primarily  of  bottom  and  deck  beams,  side  frames  and  stanchions. 
This  structure  is  called  the  transverse  ring  frame.   The  loads  this 
structure  is  subjected  to  are  side  and  bottom  water  pressure  and  the 
ship  weights  as  well  as  dynamic  loadings. 

The  structural  analysis  of  this  type  of  structure,  in  two 
and  three  dimensions,  is  discussed  in  considerable  detail  in  the 
literature.   Computer  solutions  for  both  the  two  and  three  dimensional 
aspects  are  available.   These  rapid  analysis  tools  are  necessary  pre- 
requisites for  a  rational  design  synthesis. 

In  its  most  general  form  the  optimization  of  the  transverse 
structure  entails  the  specification  of  the  transverse  frame  spacing, 
the  effect  of  longitudinal  frame  spacing,  and  the  distribution  of 
material  throughout  the  frame.   This  thesis  deals  with  the  last  part 
of  the  process,  the  optimium  distribution  of  material  throughout  the 
ring  frame  of  a  wall  sided  ship  once  the  transverse  and  longitudinal 
frame  spacings  are  fixed.   The  optimum  structure  is  defined  as  the 
least  weight  combination  of  deck  beam,  side  frames,  etc.,  that  will 
carry  the  design  loads. 
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II.   MINIMUM  WEIGHT  STRUCTURES 

The  end  and  field  moments  that  a  loaded  bean  element  develops 
are  directly  related  to  its  end  restraint.   These  moment  distributions 
are  subject  to  wide  variations.   The  bending  moment  diagrams  for  a 
beam  element  supporting  a  distributed  load  are  reproduced  in  Figure  1, 

Full  rotational  restraint  at  both  ends  produces  a  symmetri- 
cal moment  distribution  and  one  in  which  the  maximum  monents  occur 
at  the  ends,   The  field  moment  dominates  when  the  moment  restraint 
at  the  ends  is  zero.   The  intermediate  restraint  case  is  character- 
ized by  a  symmetrical  moment  distribution  and  identical  end  and 
field  moments.   In  every  case,  limiting  or  design  stresses  are 
reached  only  at  discrete  points  in  the  element. 

The  size  or  weight  of  any  element  depends  on  the  maximum 
moment  developed  when  the  design  load  is  applied,   The  wide  variation 
in  maximum  moments  is  accompanied  by  corresponding  variations  In  the 
weight  per  unit  length  of  members.   Table  1  lists,  for  the  ideal 
section  developed  in  Appendix  C,  the  variation  in  weight  for  the 
structural  configurations  of  Figure  1,   The  absolute  minimum  weight 
per  unit  length  occurs  when  the  maximum  stress  occurs  at  three  points 
simultaneously  in  the  beam  element,   The  specific  weight  increases 
by  thirty  percent  going  from  the  fixed-fixed  beam  to  the  free-free 
or  fixed-free  beam.   Doubling  the  maximum  moment  for  this  ideal 
section  increases  the  section  weight  by  sixty  percent. 
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The  weight  ratios  of  Table  1  underline  the  fact  that  the 
total  weight  of  any  structural  system  is  related  to  how  the  individ- 
ual elements  interact,   A  minimum  total  weight  of  structure  can  be 
obtained  by  substituting  members  whose  moduli  satisfy  the  inter- 
mediate restraint  condition.   In  nest  cases,  the  actual  moment 
distributions  that  obtain  when  the  structure  is  assembled  will  not 
be  identical  to  those  assumed  conditions.   The  final  total  structural 
weight  will  depend  on  how  the  material  distribution  has  been  altered 
to  obtain  convergence.   This,  in  turn,  depends  on  the  design  criteria. 
The  examples  that  follow  illustrate  how  the  total  weight  varies  when 
the  design  criteria  is  that  design  stress  must  exist  in  at  least  one 
element;  the  stress  in  all  other  members  must  be  equal  to  or  less  than 
the  design  stress. 

The  first  example,  Figure  2,  contains  only  elements  sub- 
jected to  bending.   The  length  to  depth  ratio  of  each  element  is 
assumed  to  be  large  enough  to  exclude  shear  effects. 

An  initial  choice  of  members  based  upon  the  intermediate. 
restraint  condition  and  whether  the  elements  carry  an  external  load 
gives: 

1/16  W  L2 

Z2  "    an  '    J   Zl  =  °  =  Z3 
Design 

The  structure  made  up  of  these  members  does  not  possess  adequate 

strength.   The  maximum  stress  in  member  two  is  greater  than  the 

design  stress.   The  minimum  weight  members  do  not  form  an  acceptable 

structure. 
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This  structure  can  only  be  brought  up  to  an  acceptable 
level  of  strength  by  adding  material,  The  weight  of  the  final 
structure  satisfying  the  design  criteria  will  vary  considerably  de- 
pending on  where  the  material  is  added,   There  are  two  directions 
available.   One  is  to  add  material  so  that  the  moments  at  the  nodes 
are  very  small.   This  will  cause  the  field  moment  in  member  two  to 
predominate.   The  other  members  will  be  very  smallt   The  second 
approach  is  to  increase  the  size  of  all  the  members.   This  will 
bring  about  redundant  moments  at  the  nodes  which  are  substantial, 

The  weight  curves  of  Figure  3  illustrate  the  trade  off 
possibilities  that  exist  when  all  the  members  are  varied,   Large 
increases  in  the  specific  weight  of  members  one  and  three  are 
accompanied  by  relatively  small  increases  in  member  one  as  the  dis- 
tribution factor  approaches  one.   These  effects  are  summarized  in 
the  total  weight  curve.   These  curves,  for  all  the  cases  looked  at, 
have  definite  mini -urns.   When  L  =  5L  ■  5L  ,  the  curve  shown  in 
Figure  3,  the  minimum  total  weight  occurred  when  all  member  sizes 
were  equal  and  the  node  moments  were  maximum.  At  this  minimum  point, 
design  stresses  were  obtained  in  all  members. 

The  following  example  is  identical  to  the  preceding  except 
that  all  members  have  to  support  an  external  load.   The  loading  is 
such  that  the  fixed  end  moment  of  the  middle  span  is  twice  that  of 
the  adjacent  members;  load  symmetry  is  maintained. 
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As  before,  the  structure  made  up  of  least  weight  members, 
corresponding  to  the  intermediate  restraint  condition,  is  not 
satisfactory.   All  members  are  cverstressed ,   Sufficient  strength 
can  only  be  obtained  by  adding  material. 

Total  weight  curves  for  this  structure  are  shown  in 
Figure  4,   Each  of  these  curves  have  minimum  points  which  occur  at 
k  =  .5.   At  this  point,  all  members  are  of  equal  size  and  the  node 
moments  govern.   It  is  obvious  then  that  the  minimum  total  weight 
occurs  when  design  stress  levels  exist  in  all  the  members.   At  all 
other  values  for  the  distribution  factor,  one  or  two  members  have 
too  much  material. 

Only  bending  stresses  were  considered  in  these  examples, 
The  substitution  of  a  vertical  member  at  each  of  the  supports  re- 
quires that  direct  stresses  be  accounted  for  (see  Appendix  B) ,   The 
total  weight  curve  for  this  structure  also  has  a  unique  minimum. 
This  pinimum  occurs  when  design  stresses  exist  in  all  members. 

The  iteration  technique  developed  in  this  thesis  conver- 
ges on  a  structure  which,  when  the  design  load  is  applied,  develops 
design  stresses  in  each  member  simultaneously.   These  stresses  may 
develop  at  the  ends  of  the  members  or  at  any  point  along  the 
members.   The  structures  which  this  design  process  converges  to 
represents  the  optimum  distribution  of  material  for  the  given  loads 
and  the  fixed  frame  spacing.   The  total  weight  of  these  structures 
is  a  minimum.   Any  other  combination  of  members  will  result  in  a 
higher  total  weight  of  structure;  the  stress  in  some  or  all  the 
members  will  be  less  than  tl 
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TABLE  1 


Maximum  Moc 

tent 

Position 
End,  Field 

Mome nt  Ratio 

H 
o 

Weight  Ratio 

1/16  WL2 

W 

o 

1/12  WL2 

End 

1.33  M 
o 

1.21  W 
o 

1/8  KL2 

End  or  Field 

2.0  M 

1.59  W 

WEIGHT  AND  MOMENT  RELATIONS  FOR  THE 
BEAM  SUPPORTING  A  DISTRIBUTED  LOAD 
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STRUCTURE  FOR  EXAMPLE  2,  SECTION  II 


BEAM  SYSTEMS  FOR  EXAMPLES  ONE  AND  TWO 
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III.   CONVERGENCE  PROCEDURES 

The  properties  of  a  two  dimensional  transverse  ring 
frame  of  a  wall  sided  ship  may  be  considered  similar  to  those  of  a 
plane  frame.   The  members  of  this  plane  frame  are  subjected  to  bend- 
ing and  direct  stresses;  side  sway  is  not  considered.   In  its  most 
complex  form,  this  space  frame  is  made  up  of  many  bays  and  decks. 

Picking  member  sizes  that  are  adequate,  i.e.  possess 
enough  strength  to  support  the  applied  loads,  can  be  accomplished  in 
a  straight  forward  manner,   St,  Denis  [1],  presents  a  rational 
approach  to  this  problem.   The  refinement  of  these  initial* choices, 
i.e.  cutting  the  strength  of  each  member  until  it  just  meets  the 
design  stress,  is  discussed  in  this  section, 

A  two  dimensional  space  frame  of  more  than  one  deck  is  a 
statically  indeterninent  structure.   The  stress  distribution  result- 
ing from  the  design  loads  depend  upon  the  size  of  each  member  of  the 
frame.   Changes  in  these  stress  distributions  are  expected  if  the 
size  of  one  member  is  altered.   The  effects  of  changing  the  size  of 
a  member  in  a  reasonably  balanced  structure  will  be  greatest  in  the 
immediate  vicinity  of  the  member.   It  is  assumed  in  the  following 
that  the  stress  or  moment  distribution  in  a  member  is  affected  only 
by  changes  in  members  that  join  the  defining  nodes  of  the  member. 
This  is  the  same  as  saying  that  each  member  of  the  frame  is  the 
same  as  member  2  of  Figure  5,   The  stress  distribution  in  member  2 
is  affected  by  size  changes  in  any  of  the  members  shown, 
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The  refinement  problem  Is  this:  given  a  frame  with  a  fixed 
layout  and  subjected  to  specified  design  loads,  what  distribution 
of  material  will  produce  stress  distributions  such  that  in  each 
member  the  maximum  stress  equals  the  design  stress.   This  means 
that  in  each  member  that  has  a  maximum  stress  less  than  design  stress, 
the  section  is  too  heavy.   Material  must  be  taken  away  or  the  stress 
distribution  altered  so  that  the  member  carries  more  of  the  load. 
If  the  maximum  stress  in  the  member  is  greater  than  design  stress, 
the  member  is  undersized  and  the  opposite  holds. 

An  approach  to  the  problem  that  is  suitable  to  machine 
applications  is  a  selective  procedure.  After  each  stress  analysis, 
the  member  in  the  worst  condition,  over  or  under  design  stress, 
is  selected,   The  size  of  the  member  is  altered  based  upon  existing 
maximum  moment.   The  structure  with  the  corrected  member  undergoes 
another  stress  analyses.   This  cycle  is  repeated  until  all  members 
are  at  the  design  stress  level, 

A  more  appealing  approach  is  one  in  which  changes  can  be 
made  in  all  the  members  before  each  stress  analysis.   The  results 
of  each  stress  analysis  are  used  as  the  basis  for  changing  unsatis- 
factory members. 

The  implementation  of  this  series  approach  requires  that 
.interaction  effects  be  estimated.   Unless  account  is  taken  of  how 
the  stress  distribution  changes  when  the  distribution  of  material  is 
altered,  the  system  will  be  subjected  to  undesirable  oscillations. 
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The  member  shown  in  Figure  6  depicts  a  general  member  taken 
from  the  two  dimensional  plane  frame.   It  carries  a  distributed 
load,  W,  and  has  non-symmetrical  end  restraint.   The  stress  dis- 
tribution is  not  symmetrical;  node  one  is  assured  to  be  overstressed. 

The  information  required  is  how  will  this  moment  distribution 
change  when  any  or  all  of  the  members  involved  change.   At  each  node 
three  conditions  have  to  be  considered,  an  increase,  decrease,  or  no 
change  in  the  joint  rigidity. 

An  increase  in  the  size  of  all  the  members  adjacent  to  node 
two  will  provide  additional  rotational  restraint.   This  is  accompan- 
ied by  an  increase  in  the  moment  at  node  two  ax^.d  a  decrease  in  the 
maximum  moment.   If  the  changes  in  the  members  adjacent  to  node  two 
were  sufficiently  large,  no  change  in  member  one  would  be  required 
due  to  the  altered  stress  distribution. 

Increasing  the  size  of  member  one  alone  will  cause  the 
maximum  moment  to  increase.   If  the  size  of  the  new  member  is  picked 
based  on  the  existing  maximum  moment,  then  it  also  will  be  over- 
stressed  when  the  structure  is  reassembled. 

Table  2  is  a  summary  of  the  conditions  that  have  to  be 
accounted  for.   Column  three  has  the  information  needed  to  make,  a 
rational  choice  of  new  members.   An  expected  increase  In  the  maximum 
moment  indicates  that  a  larger  than  usual  correction  should  be  made 
in  the  already  overstressed  member.   This  would  be  accomplished  by 
adding  an  increment  to  the  existing  moment,-  The  new  member  size  would 
reflect  this  added  increment, 
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Additional  tables  would  be  required  for  each  of  the 
following:  a  field  moment  governing,  end  moment  identical  and  govern- 
ing, and  non  loaded  elements. 

The  Hardy  Cross  method  of  ronant  distribution  offers  a 
direct  method  for  obtaining  numerical  estimates  of  changes  in  maxi- 
mum moments  due  to  shifts  in  material  distribution  about  nodes.   The 
expression  that  is  derived  below  actually  embodies  the  tables  dis- 
cussed in  the  last  section, 

A  first  approximation  to  the  moment  il   of  Figure  5  can 
be. obtained  by  carrying  out  the  Hardy  Cross  method  of  moment  dis- 
tribution through  the  first  unlocking  of  all  eight  joints.   It  is 
not  necessary  to  require  that  the  joints  three  through  eight  be 
fixed. 

First  release  the  joints  adjacent  to  node  two  and 
distribute  the  unbalanced  moments,  U,   The  carry  over  from  joints 
four,  six  and  eight  to  node  two  is: 

CO   -  -  1/2  k52  U6  -  1/2  k82  Ug  -  1/2  kw  U^ 

This  carry  over  is  added  to  the  unbalanced  moment  at  node  two,  U_, 
The  amount  of  this  summation  that  carries  over  to  node  one  after 
node  two  is  unlocked  is: 

1/2  kn  [1/2  (k62  U(.  +  k8,  U8  +  k^  u4>  -U2] 

The  carry  over  from  the  joints  adjacent  to  node  one  is 
found  in  the  same  way.   It  is; 

CO  -  -  1/2  k51  U5  -  1/2  k?1  U?  -  1/2  k31  U3 
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The  moment  of  node  one  in  member  1-2,  M   ,  can  be 
expressed  in  the  following  approximate  form: 

M12  '     rE''l-2  "  k12  Ul  +  1/2  k12  (k51  °5  +  k31  U3  +  k71  V 


1/2  k21  U2  ♦  1/2  k12  kH  U2  +  1/4  k21  (k62  U6  + 

k82  U8  +  k42  V 


1/4  k21  k12  (k62  u6  +  k82  u8  +  k42  u4)  +  i/4  k12  k21  u1 

1/4  k122  k21  J, 


(1) 

If  H   is  the  maximum  moment  in  the  member,  the  change  in 

the  moment  as  the  adjoining  members  and  member  two  itself  changes  is 

desired.   The  quantities  FEM.  .  and  U .  depend  only  on  the  loads  and 
n  ij      i   r 

the  geometry  and  therefore  are  constants.   The  expression  for  M 
can  be  looked  upon  as  a  function  of  the  distribution  factors.  k..» 
which  in  turn  depend  directly  on  the  section  properties. 


H12   "   f  (k12  V 


The  change  in  M   is 


dMi2  -  w~   dki2 +  ir:   dkij 

12  ij 


or 


AM.,      =     -™       Ak._    +     ~~f~       Ak., 

12  3k  12  dKi,  ij 
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Applying  this  to  equation  (1)  results  in: 


AM12  •  -  Bl  4k12  +  1/2  (k51  U5  +  k31  U3  +  k?1  0  >  ik12 


+  1/2  k^  (4k51  U5  +  Ak31  U3  +  4kn  B7)  -  1/2  ^   Ak 


+  1/2  D2  (k12  Ak21  +  Ak12  k21)  +  1/4  (k^  U6  ♦  k82  C  + 


k42  V  4k21 


+  1/4  k21  (Ak62  U6  +  Akg2  U8  ♦  Ak42  U4)  +  1/4  UA 

(k12  Ak21  +  Ak12  k21) 

(2) 

An  adequate  comparison  of  the  orders  of  magnitude  of  the 
different  parts  of  this  expression  requires  knowledge  of  the  joint 
unbalance.  If  it  is  assurr.ed,  as  it  was  previously,  that  the  loads 
are  reasonably  balanced  throughout  the  structure,  the  structure  of 
Figure  5  is  again  applicable.  Applied  to  this  structure  equation 
(2)  reduces  to: 

AM12  *  "  Ul  Ak12  "  1/2  U2  Ak21  +  1/2  U2  (k12  Ak2i  +  Aki2  k21} 

+  1/4  U1  (k12  Ak^  +  Ak12  k21) 

(3) 
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This  is  the  sane  as  neglecting  the  unbalance  at  joints 
three  through  eight,   The  inclusion  of  these  terms  in  a  computer 
program  is  straightforward,  but  of  course,  would  require  core 
computational  time.   This  final  approximateion  accounts  for  changes 
in  any  of  the  seven  members  involved,   A  more  convenient  form  for 
equation  (3)  is: 


AM..  ■  [-U.  (1  -  1/4  k,.)  +  1/2  k.  .  U.J  Ak.  . 
tj   l  i         ji'       ji  j J   ij 


+  1/2  [-  U.  (1  -  k, .)  +  1/2  U.  k. .]  Ak.. 


(A) 


Equation  (4)  gives  a  numerical  value  for  the  expected 


change  in  M   if  the  adjoining  members  and  member  itself  change  size. 
This  is  an  approximation  as  to  how  the  moment  distribution  in  member 
ij  is  affected  when  any  or  all  the  members  change  size.   This  is  a 
valid  approximation  for  both  loaded  and  unloaded  elements.   There 
is  no  restriction  to  members  undergoing  only  bending  stress;  it  is 
equally  useful  for  members  with  combination  stresses.   If  different 
types  of  materials  are  used,  their  effects  are  included  in  the 
distribution  factors. 

The  use  of  this  expression  to  find  changes  in  the  field 
moment  is  a  little  more  involved.   Equation  (4)  has  to  be  evaluated 
for  each  end  of  the  member  in  question  if  the  moment  distribution 
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is  not  symmetrical  and  added  to  the  existing  end  moments.   These 
new  moments  combined  with  the  existing  end  shears  give  an  approxi- 
mation to  the  change  in  field  moment, 

A  word  about  signs  is  appropriate.   Clockwise  moments 
are  positive t    counterclockwise  moments  are  taken  to  be  negative. 
The  signs  that  result  when  AM. .  is  calculated  cannot  be  interpreted 
literally,   A  positive  value  for  AM. .  may  indicate  an  expected  de-' 
crease  in  the  moment,  M , . ,   It  is  necessary  to  check  the  sign  of 
the  existing  moment  and  then  interpret  the  sign  of  AM  . . 
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FIGURE  6 


-  20  - 


TABLE  2 


Cumulative  Effects  of  Member  Variations 


End 

w/ 

Mom)  . 
Mxn 


End 

w/ 

Mom) 


Max 


(2) 

(1) 

Member 

(3) 

Adjacent 

to  be 

Exp,  Change 

Members 

Changed 

in  Max.  Mom, 

- 

•  + 

Increase 

+ 

+ 

Nc  Change 

o 

+ 

Increase 

- 

+ 

Decrease 

+ 

+ 

No  Change 

o 

+ 

Decrease 

Explanation  of  Table 

This  table  is  drawn  up  for  Member  Two  of  Figure  Five,   It  is 
assumed  that  the  moment  distribution  is  nonsymmetrical  and  that 
Member  Two  is  overs  tressed,   Column  three  gives  the  expected  change 
in  the  maximum  moment  based  upon  the  cumulative  effects  of  member 
changes  summarized  in  columns  one  and  two. 

Column  1:   The  sign  indicates  how  the  overall  stiffness  of 
the  joint  is  expected  to  change  based  upon  the 
stress  conditions  in  the  adjacent  members. 
Column  2:   The  size  of  Member  Two  is  expected  to  increase 
Column  3:   The  expected  change  in  the  maximum  moment  based 
upon  the  data  of  columns  one  and  two. 
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IV.   THE  DESIGN  SPIRAL 

The  result  of  the  previous  section  is  a  technique  for 
estimating  moment  changes.   These  moment  changes  are  brought  about 
when  the  members  at  a  node  change  their  size. 

This  is  of  use  in  the  design  process  because  it  presents 
a  method  of  estimating  not  only  the  magnitude  of  changes  in  the 
stress  distribution  (maximum  moment) ,  but  also  the  direction  of 
the  changes.   It  permits  a  more  rapid  convergence  to  the  final 
structure  because  each  member  of  the  structure  can  be  modified  (if 
necessary)  on  each  pass.   And  also  because  the  members  are  picked 
based  upon  conditions  that  './ill  prevail;  not  those  of  the  last 
stress  analysis.   It  is  essentially  a  method  of  predicting  the 
final  moment  distribution  and  picking  the  members  based  on  this 
distribution.   If  the  moment  changed  predicted  were  exact  instead 
of  approximations  there  wouldn't  be  any  iterative  procedure 
necessary, 

A  desirable  aspect  of  this  technique  is  that  it  is  self 
modifying  when  applied  repeatedly.   The  predicted  moment  changes 
are  updated  each  time  a  new  member  is  picked.   In  this  way  oscilla- 
tions are  avoided  and  smooth  convergence  is  assured.   It  is  not 
necessary  to  perform  a  new  structural  .analysis  until  this  self 
modifying  aspect  is  completed.   This  occurs  when  there  are  no 
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significant- differences  in  the  predicted  moment  changes  from  one 
cycle  to  the  next. 

The  rational  design  process  developed  in  the  sections  that 
follow  converges  to  the  opt     distribution  of  material  for  a  trans- 
verse ring  frame.   The  general  procedure  followed  is  to  use  the 
results  of  a  stress  analysis  as  the  basis  for  modifying  the  material 
distribution  in  such  a  manner  that  design  stresses  will  be  attained 
in  each  member,   The  results  of  Section  II  indicate  that  a  structure 
which  develops  design  stresses  in  each  member  simultaneously  is 
the  least  weight  structure. 

Equation  (4)  is  used  in  this  process  to  make  corrections 
to  the  results  of  the  stress  analysis,   A  change  in  the  scantlings 
of  any  member  of  the  ring  frame  changes  the  stress  distributions. 
As  pointed  out  in  Section  III  and  summarized  in  Table  II,  the 
effect  of  a  change  in  scantlings  depends  not  only  en  the  magnitude 
of  the  change,  but  also  the  existing  moment  distribution.   Once  the 
maximum  moment  in  the  member  is  identified,  equation  (4)  gives 
direct  information  as  to  how  this  moment  will  change  when  the 
scantlings  of  the  members  of  the  ring  frame  are  revised. 

This  correction  to  the  existing  moment  distributions 
provides  the  means  of  facilitating  rapid  convergence  to  the  optimum 
structure.   Changes  in  member  scantlings  are  based  on  these  pre- 
dicted moments.   The  objective  of  design  stress  in  each  member  of 
the  ring  frame  is  obtained  much  more  readily. 


■  The  use  of  equation  (A)  in  the  design  process  requires 
the  standard  output  of  a  stress  analysis,  end  forces,  plus  end 
moments.   The  process  begins  with  the  results  of  the  stress  analysis 
which  used  the  initial  choices  of  members. 

The  computer  begins  with  member  one,   The  governing 
moment  is  ascertained  and  the  stress  level  is  computed  based  on  this 
maximum  moment ■   If  the  stress  level  is  within  a  certain  prescribed 
range  of  the  design  stress,  the  member  is  not  changed.   This  pro- 
vides a  means  for  ending  the  iterations.   If  the  stress  level  is 
not  satisfactory,  a  new  member  size  is  calculated  and  stored  in  the 
matrix,  K,   This  procedure  is  repeated  for  each  member  of  -'the 
structure.   The  result  of  this  cycle  is  the  ratrix  K,   The  members 
that  make  up  K  represent  the  first  cut  at  bringing  about  the  desired 
distribution  of  material.   These  members  are  fictitious  in  the 
sense  that  they  are  not  used  in  an  actual  stress  analysis.   The 
properties  of  these  members  are  based  upon  the  existing  stress 
distributions. 

The  computer  now  starts  to  revise,  the  members  of  matrix  K; 

it  starts  with  the  first  loaded  member.   The  stiffnesses,  k..,of 

the  members  at  the  defining  nodes  of  the  -ember  to  be  changed  are 

computed  using  the  data  of  matrix  K,   This  data  is  used  to  compute 

AM , .  or  AM.,  or  both  depending  on  the  stress  distribution,   If  the 
ij      jx 

field  moment  governs  AM, .  is  computed  for  each  node;  when  an  end 

ij 

moment  governs  AM.,  is  calculated  for  the  governing  node. 
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The  maximum  moment  of  the  member  to  be  changed  is  modified 
by  adding  the  correction  factor  AM...   The  properties  of  a  new 
section  are  calculated  based  on  this  modified  moment.   These 
properties  are  inserted  in  the  proper  slot  of  matrix  K, 

The  computer  chooses  the  next  loaded  member  and  repeats 
the  above.  The  unloaded  members  are  revised  after  new  properties 
have  been  picked  for  all  the  loaded  members.  This  requires  fewer 
steps  because  the  end  moments  always  govern. 

The  result  after  all  the  members  have  been  looked  at  is 

i 

a  new-  matrix  of  members,  K  ,  A  stress  analysis  could  be  initiated 

at  this  point.   The  results  of  the  examples  indicate,  however,  that 

a  smoother  convergence  is  obtained  if  this  process  is  repeated,  i.e. 

ii 

find  the  matrix  K   , 

A  stress  analysis  is  now  undertaken  using  the  members  of 

ii 

the  matrix  K   ,   The  results  of  this  stress  analysis  are  used  as 

the  starting  point  for  a  new  design  cycle. 

Flow  charts  containing  the  details  of  the  design  spiral 
are  contained  in  Appendix  A,   Illustrations  of  the  design  spiral 
in  action  are  contained  in  Appendix  B, 
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V.   RESULTS  ALT)  CONCLUSIONS 

A  rational  design  procedure  has  been  developed  which 
permits  the  rapid  synthesis  of  the  transverse  ring  frame  of  a 
wall  sided  ship.   The  resulting  transverse  structure  represents 
the  optimum  distribution  of  material,  the  minimum  wieght  combus- 
tion of  deck  beans,  side  frames,  etc,,  within  the  constraints  of 
fixed  longitudinal  and  transverse  frame  spacing. 

Detailed  applications  of  this  procedure  are  carried  out 
in  Appendix  B,   Each  step  follows  the  subroutines  developed  fcr 
computer  adaptation  of  this  technique.   Mo  attempt  has  been  made 
to  adapt  the  procedure  to  hand  calculations. 

This  method  converges  rapidly.   Large  member  changes  were 
required  in  each  of  the  examples  to  attain  the  design  stress  in  each 
member.   Example  I  required  only  one  additional  stress  analysis 
after  the  initial  choice  of  members  were  analyzed;  Example  II  re- 
quired two  additional  moment  distributions. 

The  speed  with  which  a  satisfactory  ring  frame  is  obtained 
depends  directly  on  the  following:  the  limits  set  within  which  the 
member  stress  must  fall  and  the  use  of  weighting  functions  in  con- 
junction with  the  predicted  moment  changes. 

The  limits  set  for  member  stresses  in  Appendix  B  were 
liberal.   This,  of  course,  was  conducive  to  a  rapid  convergence. 
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These  limits  are  not  fixed,  however,  and  can  be  set  at  any  desired 
level.   It  is  possible  to  put  new  limits  on  the  members  stress  with 
each  problem. 

The  predicted  moment  changes  in  the  examples  were  used 
as  given  by  equation  (4) ,   The  true  approximation  in  equation  (4)  is 
in  the  magnitude  of  expected  moment  change.   The  direction  of  the 
change,  plus  or  minus,  is  exact  if  the  change  is  greater  than  a  few 
percent.   In  all  the  cases  investigated,  the  magnitude  of  the  change 
was  always  low.   Convergence  would  have  been  more  rapid  if  the 
changes  were  multiplied  by  a  weighting  factor.   This  indicates  that 
weighting  factors,  if  necessary,  can  be  developed  as  experience  is 
gained  with  the  basic  program.   These  weighting  factors  would  become 
part  of  the  basic  program. 

The  results  of  Example  II  of  Appendix  B  indicates  that 
members  whose  continuity  is  not  required  to  carry  the  applied  loads 
may  be  eliminated.   This  is  easily  avoided  by  applying  nominal 
loads  to  these  members,  e.g.  their  own  weight.   In  an  actual 
structure,  this  will  seldom  be  a  problem. 

The  use  of  the  ideal  section  of  Appendix  C  was  a  necessity 
in  the  development  stages  of  this  program.   The  implementation  of 
this  design  procedure,  i.e.  the  actual  programming,  does  not  re- 
quire this  ideal  section,   If  a  tape  of  actual  sections  is  available, 
the  program  can  use  it  and  still  work,   The  effect  this  would  have 
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on  the  speed  of  convergence  would  be  mixed.   The  actual  members  in 
most  cases  have  larger  section  moduli  than  those  the  program  would 
pick,  "  The  use  of  actual  sections  requires^in  addition,  the  develop- 
ment of  a  trial  and  error  solution,  in  the  subroutine  COMSTRESS, 


r 
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APPENDIX  A 


An  algorithm  covering  the  complete  design  process  is 
developed.   Flow  charts  covering  the  details  of  the  individual  steps 
of  the  algorithm  are  also  presented. 

Design  Algorithm 


Start  with  member  one  compare  the  maximum 
and  design  stresses.   If  the  maximum  stress  is 
within  the  prescribed  range,  insert  the  pre- 
sent member  properties  in  matrix  K,   If  maxi- 
mum stress  unsat,,  pick  a  new  member,   Insert 
its  properties  in  the  matrix  K, 


Repeat  the  above,  for  each  member.   If  all 
members  have  satisfactory  stress  levels  the  com- 
puter stops  -  the  design  is  completed.   If  any 

or  all  members  were  changed,  the  iteration 
process  is  begun  by  returning  to  the  loaded 
member  number  one. 


Required  stiffness  factors  are  calculated 
using  elements  in  matrix  K,   The  moment  change 
is  determined  using  this  data.   The  existing 
maximum  moment  plus  this  moment  change  is  the 
new  I     .  '  :   nt,  k   new  :  - 

this  design  moment  and  the  direct  stress. 
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The  new  member  is  placed  in  matrix  R, 
This  process  is  repeated  for  all  members 
even  if  stress  levels  are  initially  satis- 
factory.  The  matrix  containing  these  new 

i 
members  is  now  designated  K  , 


The  data  of  matrix  K  is  now  used  to 

calculate  the  required  stiffness  factors. 

The  process  is  repeated.   The  matrix 
ii 

formed  is  K   , 


A  stress  analysis  is  carried  out  using 
matrix  K   ,   The  results  of  this  analysis 
are  the  starting  point  for  a  new  design 

eye'  : , 


The  flow  charts  are  made  up  in  the  form  of  subroutines. 
They  cover  all  the  steps  presented  in  the  design  algorithm. 

The  first  flow  chart  is  subroutine  Moment.   It  uses  as  an 
input  the  data  of  the  preceding  stress  analysis.   This  subroutine 
primarily  finds  the  maximum  moment  and  calculates  the  maximum 
stress.   When  field  moments  govern  the  member  is  given  a  special  tag. 

Subroutine  Change  checks  the  maximum  stress  in  each  member. 
It  stops  the  iteration  process  if  all  members  are  satisfactory. 
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Subroutine  Stiffness  uses  the  data  of  matrix  K  and  the 
data  of  initial  member  choices  to  calculate  the  expected  changes  in 
stiffness.   It  also  determines  the  moment  change  and  adds  this  to 
the  maximum  moment.   The  new  member  parameters  are  obtained  by 
calling  the  appropriate  subroutine,  i.e.  Pick  or  Comstress,   This 
subroutine  assumes  that  a  member  numbering  system  is  available  which 
facilitates  picking  the  required  data  from  the  matrix  K,  etc. 

Subroutines  Pick  and  Comstress  both  choose  new  member 
properties,   Comstress  handles  members  that  have  significant  direct 
stress  as  well  as  bending  stresses. 

The  calling  sequence  for  these,  subroutines  would  be: 

CALL  SUBROUTINE  MOMENT 

CALL  SUBROUTINE  CHANGE 

CALL  SUBROUTINE  STIFFNESS 

CALL  STRESS  -  or  any  suitable  stress  analysis  program. 
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SUBROUTINE  STIFFNESS 
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APPENDIX  B 
EXAMPLES 

EXAMPLE  I 

The  seven  membered  frame  shown  on  page  45  is  used  to 
illustrate  the  optimization  technique  developed  in  the  body  of 
the  thesis. 

A.  Selection  of  Initial  Member.  Jl£2££E£iisSi 

'  The  initial  set  of  member  properties  is  part  of  the  data 
input  to  the  computer.   The  following  were  selected  for  this 
example* 

(a)  Z  =  240  in3  ;  I  =  2000  in4 

(b)  Z±   =  Z3  =  120  in3;  I   =  I   =  320  in4 

(c)  Z.  -  Zc  =  Z^  -  Z,  -  120  in3;  I,  =  820  in4 

4    5     6     7  4 

These  properties  were  obtained  using  the  ideal  section  of  Appendix  C, 

B.  Stres     ilysis 

The  Hardy  Cross  method  of  moment  distribution  was  used  to 
obtain  end  moments  and  reactions,  The  results  of  this  first  analy- 
sis are: 

M.„  =  M_.  =  M.)M   =  684  KIP  FT 
12    21    1  Max 

M.7  =  VL     m   m  )    .  495  KIP  FT 
17    28    2  Max 

Mn  =  MA>Mav  =  95  K»;  F,  =  130  KIP  FT 


13    4 'Max 
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Cjt|  Subroutine  Change 

(1)  The  maximum  stress  In  each  member  is  compared  to 
the  design  stress  =  40  KSI 

(2)  The  size  of  each  unsatisfactory  member  is  changed. 
The  maximum  stress  in  each  member  is: 

o,)v   -  34  KSI 

1  Max 

o,)M   =  50  KSI 
7  Max 

a  )Mm  «"  15  KSI  +  10  KSI  =  25  KSI  -  a  combined  stress  member. 
The  new  properties  for  each  member  are: 

M  ) 

z     =  f-   ^,x  «   204  in  ;  I_  =  1600  in 

2  0-..  2 

Design 

Z  =  148  in3;  I  =  1100  in4 

Z.  -  67  in3;  A.  -  3.3  in2;  I.  =  415  in4 
4  4  4 

The  remaining  member  sizes  follow  by  symmetry.   These  proper- 
ties make  up  the  matrix  K,   The  subroutine  PICK  would  be  uied 
to  obtain  the  properties  of  members  one  through  three;  sub- 
routine COMSTRESS  would  be  called  to  obtain  the  properties 
for  members  four  through  eight, 

D .   Sub  rout i/rv e^  S  t_i f  f n_ e  s s 

i 

(1)   Calculates  new  stiffness  factors,  k.. 

■ '  ij 

,(2)   Calculates  the  charge  in  stiffness  factors  with 

respect  to  the  initial  member  choices. 
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D.   Subroutine  Stiffness   (Cont'd) 

(3)  Determines  AM   and  the  new  factitious  maximum  moment 

(4)  Picks  new  members  based  one  the  new  maximum  moment. 

This  subroutine  begins  with  the  first  loaded  member,  member  1 

i 

and  calculates  k17 


k1?  =   .310 

k?1  =  o 

k1?   -    .184 

kn-o 

t 
Ak17   =  +   .126 

Akn  -  0 

These  results  are  used  in  equation  (4)  to  obtain  AM-.. 


AM    =  +  50.4  K'  ;  M    -  +  495  +  50.4  =  545  KIP  FT 

Subroutine  Pick  noi;  comes  up  with 

'         3    '  4 

Z  =  153  in  j  I  =  1140  in 

These  member  properties  are  immediately  inserted  in  the 
matrix  K.   This  provides  a  constant  updating  of  the  data. 

The  computer  proceeds  thru  the  structure  doing  the 
above  for  each  member 
Member  2 

k12  =  .448       k2  -  .452 

k   =  .450       k2]  -  .450 

Ak12  =  -.002    Ak2]  =  +  .002 

AM*  2  =  1.1K';  -  M  2  =  685  KIP  FT 

Z     -  204  in3;  I  -  1600  in4 


*   JU   * 


D.   Subroutine  Stiffness   (Cont f  d ) 
Member  3 


k,,  =  .310  fc  ,  =  0 


17    "™  "71 


k17  "  a84  k71  "  ° 


:1?  -  +  .126         Ak   «  0 


i 

A 


^  =  +  50.4  KIP  FT    M_  =  545  KIP  FT 


»  T       » 


Z7  -  153  in3;  I_  ■  1140  in4 


7 

J.  -'  J         .Li. 

» 

Nenber 

4 

i 

k13  ' 

=  .117 

k13~ 

.184 

1 

4k13- 

-  .067 

t 

U±- 

27  KIP 

FT 

i 

k31  = 

•  0 

k31  = 

=  0 

1 

4k31 

-  0 

1 

S" 

73  K 

IP 

FT 

•        3    '         4 
Z  ■  55  in j  I  =  310  in      Subroutine  COMSTRESS 


Meirber  5 


k   *  .121  k51  =  0 


Ak.c  =  .184  Akcl  =  0 

15  51 


Ak15  =  -.063  k51  -  0 


AM  =  -25.2  KIP  FT     M  =  75  KIP  FT 


Z  =  55  in  ,  A  =  7,2  in  ;  I  -  310  in 


.  2   T' 

in  ; 

Members  6  and  7  follow  by  symmetry 
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» 
The  matrix  K  at  this  point  becomes  the  matrix  K  , 

Subroutine  Stiffness  repeats  the  above  process  and  produces  matrix 

n 

K  .   The  results  of  this  procedure  are  summarized  below 

Member      Z   (in3)    A   (in2)      I   (in4) 

1  167        — -  1300 

2  204 1600 

4  50       6.9  280 

5  50        6.9  280 

ii 

The  matrix  K   is  now  the  matrix  of  actual  member 

properties,   A  stress  analysis  of  the  structure  made  up  of  these 
members  is  performed ;   The  results  of  this  analysis  are: 

M.  =  595  KIP  FT     o,  =  42  KSI  =  1.07  a 

1  1  o 

M„  =  630  KIP  FT     a  =  40  KSI  =  a 

2  2  o 

M,  =  .41  KIP  FT     a,  =  36  KSI  -  ,9a 
4  4  o 

Stresses  and  moments  for  other  members  follow  by  symmetry. 

Subroutine  Change  now  finds  all  the  stress  levels 

ii 

satisfactory.   The  members  of  matrix  K  represent  the  optimum 

structure.   The  design  spiral  is  complete, 
EXAMPLE  2 

This  example  was  chosen  to  illustrate  the  procedure  for 
a  nonsymmetrical  system  and  one  in  which  the  field  moment  governs, 
at  least  initially. 
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The  frame  analyzed  is  shown  on  page  46  .   The  optimization 
process  proceeds  as  illustrated  in  detail  in  example  one, 
A,   Initial  i   ber  E  onerties 

Z2  =  240  in3;  I  =  2000  in^ 

3  2  A 

All  other  members  -  Z  ■  40  in  ;  A  =  6  in  ;  I  =  210  in 

B|ti  j.tre  ssAn  a  ly  s  i  s 

^  =  16 S  KIP  FT 
M2  =  780  KIP  FT 
M,  =163  KIP  FT 
M5  -  84  KIP  FT 
C_.  t  .SuJ^PJj.tine. .Chan ?te 

ai  «  50  KSI  Zj  =  50  in3;  I  -  290  in4 

o2  «  39  KSI  ;  My  governing        , 

Z   =  240  in  ;  I  -  200  in 
a.    =  50  KSI  L  l 

Combined  Stress    ,  ,  , 

a   =  38.5  KSI  z  -  83  in3;  A  =  9.5  in  ;  I   -  560  in' 


i        3    '         2    * 

Z5  -  37  in  ;  A  =5,6  in  ;  I  -  190  in 


D.   Subroutine  Stiffness 


1 

Matrix  K 

Member 

Z 

•    3 

(in  ) 

1 

56.1 

2 

203 

4 

140 

5 

23 

'2  '4 

A   (in  )  I   (in  ) 

320 

1700 

13.5  1050 

4.1  100 
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The  field  moment  governed  in  member  2  at  this  point. 

However  o:ice  the  second  set  of  members  were  obtained,  the  end 

moments  p redo minated.   The  results  of  the  second  pass  is 

n 
matrix  K 

"3       "2       "4 

Member     Z   (in  )     A   (in  )     I   (in  ) 

1  56,1        320 

2  203         1650 

4  190         16.5         1500 

5  15  3  60 

E.  ^  Stress  Analv5isi 


M  =  530  KIP  FT 


M2  »  655  KIP  FT 


M4  »=  114  KIP  FT 


M  =  0 
F.   Subroutine  Chanpe 


on  =  38.8  KSI  =  .97 
1  o 


a0  =  30  KSI  -  .75 
2  o 


a,  =  48  KSI  =  1.25 
4  o 


o5  =  0 
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These  stress  values  require  another  pass.   The  trend  is 
obvious.   Members  two,  four  and  six  are  taking  all  the  load. 
The  second  iteration  with  subroutine  Stiffness  produces  the 
follovinr:  matrices. 


Matrix  K  -  second  pass 


Member 

1    3 
Z   (inJ) 

'     2 
A   (inZ) 

I   (i 

1 

17 

— 

40 

2 

203 



1650 

4 

262 

21 

2000 

5 

— — 

____ 

__„„ 

Matrix  K  -  second  pass 

it     -  "2  "4 

Member     Z   (in  )  A   (in  )  1   (in  ) 

2         192  1550 

4         262  21  2000 


G,   Stress  Analysis 


M2  »646  KIP  FT 


M^  ■=  646  KIP  FT 


H.  i  Subroutine  Change 


o-  =  39.5  KSI   -   .98  o 
2  o 

a,    -  29.6  +  11.1  -  40.7  KSI  =  1.025  o 
4  o 
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These  stresses  are  within  the  specified  limits, 

it 

The  design  spiral  is  terminated.   The  members  of  matrix  K  , 

second  pass,  form  the  optimum  structure. 


-  44  - 


\\\ 

\\  \ 
3 

4 

24  KIPS/FT 

© 

L  -  20» 

12  KIPS/FT 

V 

'                   v 

t 

12  KIPS /FT 

- 

^ 

f 

h       © 

2     © 

2       ©         8 1: 

© 

© 

wx 

5 

AT  S 

rrr 

6 

STRUCTURE   FOR  EXAIITLE   B-l 


FIGURE   B-l 


-   45   ■ 


\  i» 

\  \  \ 
3 

l_Y  V 

A.  V.V 

4 

© 

24  KIPS /FT 

© 

i 

\ 

' 

'i 

' 

> 

i  i 

® 
© 

1                    © 

2 

© 

VTT 

5 

T-T- 

6 

Ll  ■   L3   =  L4 


L5   »    10' 


L0   -  hr    »  U   -  20' 


STRUCTURE   FOR  EXAMPLE   B-2 


FIG'JRE   B-2 


-  46  - 


APPENDIX  C 
THE  IDEAL  SECTION 

The  section  moduli  of  available  steel  sections  listed  in 
t*ie  Manual  of  S^teel  Construction  do  not  vary  in  a  continuous 
manner.   The  use  of  the  least  weight  sections  compounds  the 
problem.   Each  section  is  useful  for  a  wide  range  of  section  moduli. 

These  facts  make  it  very  difficult,  if  one  is  using  the 
sections  listed  in  the  AISC  Manual,  to  pick  out  trends.   Also  the 
fact  that  there  are  no  mathematical  relationships  existing  between 
the  important  properties  of  the  sections  compounds  the  problems 
encountered  when  using  actual  sections, 

A  mathematical  section,  i,e,  one  in  which  there  are 
specific  relationships  existing  between  the  section  parameters, 
avoids  these  problems.   The  section  developed  for  this  thesis  is 
outlined  below.   It  conforms  well  to  minimum  weight  sections  listed 
in  the  AISC  Manual  for  larger  section  moduli.   At  the  lower  end  of 
the  scale  there  is  a  greater  divergence  from  actual  sections, 

Let  tf  ■  flange  thickness 

t  ■  web  thickness 
w 

D  =  depth  of  section 
wf  ■  flange  width 

The  relationships  chosen  for  the  section  are: 

t,  -  2  t 

f      w 

D  =  50  t 
w 

w.  =  20  t 
f       w 
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The  expressions  for  the  area,  section  modulus  and  moment 
of  inertia  are: 

A  =  .052  D2 
Z  =  .032  D3 
I  -  .0154  D4 

A  log-log  plot  of  these  expressions  is  made  on  page  49. 
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